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^3 ' Abstract 

a: 

Q-i! We analyze oscillations of intermediate neutrinos in terms of the scattering of 

particles described by Gaussian wave packets. We study a scalar model as in a previous 
paper (I) but in realistic situations, where the two particles of the initial state and final 
state are wave packets and neutrinos are in the intermediate state. The oscillation 

b : 

of the intermediate neutrino is found from the time evolution of the total transition 
probability between the initial state and final state. The effect of a finite lifetime and a 
finite relaxation time are also studied. We find that the oscillation pattern depends on 
the magnitude of wave packet sizes of particles in the initial state and final state and 
the lifetime of the initial particle. For Z\m|i = 3 x 10~ 2 eV 2 , the oscillation probability 
deviates from that of the standard formula if the wave packet sizes are around 10~ 13 
m for 0.4 MeV neutrino. 



§1. Introduction 

Neutrino oscillation is the only phenomenon in which to see effects of neutrino masses 
at present. To analyze neutrino oscillations, single particle wave functions have mainly been 
studied. Because neutrino masses are very important, it is necessary to understand the 
quantum mechanics of neutrino oscillations in full detail. It is the purpose of the present 
paper to study quantum mechanical aspects of neutrino oscillations beyond the single particle 
picture. We studied particle oscillations from a nonstandard viewpoint in the framework of 
quantum field theory, where neutrinos are in the intermediate state and the finite time 
interval effect is explicitly taken into account, on the basis of plane waves in a previous 
paper. ^ We found that non-standard oscillation patterns emerge in the exact plane waves. 
In many real physical processes, however, particles are not exact plane waves but have finite 
spatial extensions. Wave packets are suitable to express these particles. We study particle 
oscillations in terms of the scattering amplitude of the particles described by Gaussian wave 
packets. Neutrinos are in the intermediate state in this amplitude, and neutrino oscillation 
is studied from this amplitude. 

In scattering processes, the roles and importance of wave packets have been stressed by 
Goldberger and Watson. 2 ) The size of the wave packet of the initial state is determined from 
the beam size, and it is a semi-micro scale for hadron beams. This corresponds to an energy 
scale of the order of eV or less. By contrast, typical energy scales of hadron systems are of 
the order of a few hundred MeV. Hence the effects of wave packets in the initial state are 
negligible in standard hadron experiments. The size of the wave packet of the final state, 
on the other hand, is determined by a detector. A detector is composed of many systems 
of materials. We regard the minimun set of materials for which a classical signal is taken 
as a unit detector. A unit detector is composed of atoms and generates radiation, electrons, 
or other particles by which information from the quantum wave function is transmitted 
to classical observers. The wave packet size is determined from the unit detector and the 
resolution of the measurement is also determined from the unit detector. Hence the wave 
packet size should be about the same as the spatial resolution of the detector. A unique 
value of this size is not known now, so we study the dependence of the transition probability 
on the wave packet sizes. The time resolution is determined from the time evolution of the 
detector. In this paper, for simplicity, it is assumed that the time resolution is zero. Two 
observations at different times are assumed to be independent. 

In neutrino oscillation experiments, the typical energy scale is extremely small and the 
spatial sizes of experiments are of the order several hundreds km or more and are very 
different from scales of ordinary experiments. Hence, the roles of wave packets in neutrino 



experiments are different from their roles in other ordinary experiments and they should 
be clarified. It is one of our purposes in the present work to study these problems. In 
fact, several theoretical works have been done on Gaussian wave packets. But they are not 
sufficient. Especially, qualitative analysis is lacking. We study effects that have not been 
studied to this time. The works in Refs. 3)-5) treat neutrinos as wave packets in a single 
particle picture. Other works treat the particles in the initial and final states as wave packets 
in a field theoretical treatment. 6 )' 7 - 1 In the latter, the standard S matrix theory, in which the 
transition time interval is set to infinity from the beginning, was used, and the dependence 
of the amplitudes on neutrino parameters was obtained. However, this standard treatment 
of the S matrix is inadequate in a process in which the finite time interval effect is important 
and the dependence of amplitudes on external particles' parameters are studied. The finite 
time interval effect becomes relevant in the situation in which the intermediate particle 
is very light and interacts weakly with matter. Especially when the intermediate particle 
consists of a superposition of several mass eigenstates and the mass squared difference is 
very small, the finite time interval effect is not negligible. In the standard S matrix theory, 
because the energy is strictly conserved and the interference of the amplitudes of different 
mass does not occur, there is no oscillation when all the particles are exact plane waves. In 
(I), the finite time interval effects in particle oscillations was shown to be important when 
the observed particles are exact plane waves. Thus the finite time interval effect should be 
important generally in the field theoretical treatment of neutrino oscillations. A modified S 
matrix approach that allows us to investigate finite time interval effects should be applied 
for the study of intermediate neutrino oscillations. 

In the present paper, we extend our study of particle oscillations in the intermediate 
state to the wave packet formalism. Oscillation amplitudes of neutrinos in the intermediate 
state where the particles in the initial state and final state are described by wave packets 
are studied and effects due to the finite wave packet, the finite time interval, and the finite 
lifetime or finite relaxation time are found. Amplitudes are shown to deviate from the 
standard formula in extreme conditions, when the wave packet sizes are very small. 

This paper is organized in the following manner. In Section 2 we give the general con- 
sideration of the wave packet formalism in which the particles in the initial and final states 
are described by Gaussian wave packets in the finite time interval method. In Section 3, the 
amplitude is computed using the Gaussian approximation. In Section 4, we include a finite 
lifetime and a finite relaxation time. In Section 5, numerical results in one spatial dimension 
are presented. A summary is given in Section 6. 



§2. The wave packet formalism 

Here we investigate neutrino oscillations in a scalar model in which particles in the initial 
and final states have finite spatial widths that are described by Gaussian wave packets. 
Particles A, B, C and D are external particles and are expressed by the field operators 
<H>a(x), @b(x), @c(x) an d &d(x). The fields ^(x) and &i 2 (x) are mass eigenstates and are 
internal particles. The Lagrangian density is given by 






(2-1^ 



where m^ is the mass of <Pi i and mi that of <?£,. The interaction Hamiltonian is written 

Hl nt = F x fd 3 x <P A (x)<P Ic (x)<P c (x), 

Hf nt = F 2 fd 3 x <P b (x)<P Id (x)<P d (x), (2-2) 

where F\ and F 2 are coupling constants. The fields <Pi c (x) and <Pi D (x) in the above interac- 
tion Hamiltonian are linear combinations of the mass eigenstates : 

$i c (x) = cos 9 ■ <P h (x) + sin 9 ■ $ l2 (x) , 

&i D (x) = -sin 9 -^(x) + cos 9 -<P l2 (x). (2-3) 

Here 9 is the mixing angle between <Pj c and <Pi D and between <Pj 1 and <Pi 2 . 
Each field operator is expanded in the interaction representation as 

$l{x) = / . , expf ip ■ x - iE L (p)tj ■ a(p) L + h.c, 

E(p) L = vV + ^l 2 , (2-4) 

where L stands for A,B,C,D and ij. 

We investigate the situation in which the particles A and B are prepared at time t = 
and position x = X^, and the particles C and D are detected at t — T c , x = X c and 
t = T D , x = Xx>, respectively. (Fig. |TJ) We assume To > T C - The transition amplitude of 
the finite time interval Tjj between the initial and final states at second order is calculated 
as 

(final) S[t = T D} t = 0] (initial) = (final) i 2 / dt 2 / d£i# int (£i)# int (£ 2 ) [initial). (2-5) 
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Fig. 1. Diagram 

The initial state is composed of two particles, A and B, which are the wave packets of 
the finite spatial extents expressed by a distribution function. We assume that the initial 
time and the final time are defined with infinite precision. This assumption makes the wave 
packet states defined at different times independent of each other. The distribution function 
of the momentum p, w;(p ; p , X,T; a), has a finite extension around the central value p°, 
and a Gaussian form of width a is assumed : 

-ip-X + iE(p)T 



w(p;P°,X,T;or) = exp 
The initial state is given by 
| initial) = 
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(2-6) 
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The final state is composed of two particles, C and D, which are the wave packets. The final 
state is defined in the form 

d 3 p c f d 3 p D 
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(2-8) 



In the above equations, X^, Xg,Xc and X^ are the center positions of the wave packets, 
and p^, p'g, p^ and p^ are the central values of momentum of the particles A, B, C and 
D, respectively. We consider the spatial sizes of the wave packets, a xL = 1/(2o"l) (L = 
A, B, C, D), to be between a macroscopic size and a microscopic size. 

In the amplitude Eq. (J2-5J1 . the particles A, B, C and D represent directly observed 
particles. The particle I represent a scalar neutrino and appears only in the intermediate 
state. Substituting H int (t), the amplitude is given by 



S = (final) S[t = T D ,t = 0] (initial) = -F X F 2 sin 26 (S x - S 2 ) , 



(2-9) 



where 
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(2-10) 



From this amplitude, we study neutrino oscillation in the intermediate state. Because the 
neutrino interacts with matter extremely weakly, it is not observed directly in real exper- 
iments. Hence, the amplitudes of the present situation agree with the amplitudes of the 
realistic experimental situations. 

Now we perform the integrations over p^, p#, pc and po in Eq. (J2-10J1 . To integrate 
these variables, the energy is expanded around its central value, 



E L (p) = E L (p L ) + (p-p L )-v L , 
dE L (p° L ) 



vl 



dpi 



(2-11] 



where v^ is the velocity of the particle L at this momentum, and the momentum integration 
is carried out as Gaussian integrations. After similar integrations over x 1; x 2 and k, the 
transition amplitude is obtained as 

AP 2 



Si = N exp 
where 

iV = 



2a 2 



dt 2 I dt\ exp 
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-Zi(h, t 2 , T D ) + ihAEl + it 2 AE° 2l 



(2-12) 
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-X A • p° A - X B • p° B + (Xc • p° c - T C E° C ) + (X D ■ p° D - T D E° D ). 



(2-13) 
(2-14) 



The quantity Zi in the exponent causes the amplitudes Si to have dominant contributions 
from the regions Zi ~ and is called a "trajectory function". The trajectory functions 
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Fi{t u t 2 , T D ) = x°(£ 2 , T D ) - x?(tx) - Vi (t 2 - £i), 
G(t 2 , T D ) = X D - X B - t 2 \ B - (T D - t 2 )v D , 
H(ti) = X C -X A - £ lVA - (T c - ti)v c 



(2-15) 

(2-16) 
(2-17) 
(2-18) 



give a classical particle picture, where Vj is the velocity of Jj and is defined in the same 
way as 1)2-11)) . Here we use for conciseness a 2 = o\ + a 2 B + a c + a 2 D and a 2 - = erf + a 2 



(i, j = A,B,C,D). The coefficient of F 2 (£i,£ 2 , ^d), °\c (J 'bdI o2 i ls the width of the Gaus- 
sian function for the intermediate particle's momentum. In our approach, the intermediate 
particles become wave packets automatically and then sizes are given by those of the external 
particles. 

The momenta and energies appearing in <fi are given by, 

zAP = p° c + p° D - p^ - p% (2-19) 

a 2 
pO =p o + _Mp f or L = A,B (2-20) 



r 2 



pO=pO-^|Z\P for L = C,D, (2-21) 
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E° L = E° L + -^v L -AP, for L = A,B (2-22) 

2 

E° L = E° L -^v L -AP, for L = C,D (2-23) 

AE° U = E C + E t (k )-E° A , (2-24) 

AEl = E° D - ^(k°) - E° B . (2-25) 

The quantity — Z\P 2 /(2cr 2 ) in the exponent of Eq. (J2-12J1 is called the "momentum function". 
It gives a constraint on the differences between momenta. 

The intermediate particles appear as wave packets even if they were not originally pre- 
pared as wave packets, and their momenta are given by the following function of the momenta 
of external particles : 



P';> - pL - 4^P. (2-2G) 



x x and x° are the central positions of the interaction vertices and written 

x$(*i) = 4-{°a( X a + tiv A ) + 4(X C - (T c - ti)v c ) 1, (2-27) 

x°(t 2 , T D ) = J-L|(X B + hv B ) + 4(X D - (T D - t 2 )w D )\. (2-28) 

a BD I J 

The Gaussian term in Eq. (J2T2J1 places the constraint p^ + p° D — p A — p^ ~ on the 
momentum with the width \j2a 2 . The momentum is approximately conserved, because the 
initial state and final state are approximate eigenstates of the momentum. The Gaussian 
terms in the integrand of Si place the constraints F, m 0, G ~ and H m on to the times 
t\ and ^2- These correspond to the classical trajectories in the particle picture. The latter 
constraints become stronger as the spatial sizes of the wave packets become smaller. 



§3. Gaussian approximation in time integration 



In this section we study the amplitude Eq. (|2T2|) further. In order to perform the time 
integrations, the exponent of the integrand in Eq. (|2T2|) is rewritten as follows 

- lZi{t u t 2 , T D ) + ihAEl + it 2 AEl = ~Z t (t° u , t° 2i , T D ) + tt^T^AE , + tt° 2l (T D )AE° 2l 
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- \ (*t u AEtf - 1 (a t2i AE° 2i + At\ t AEl)\ 
where a\ u and af 2 . are the widths of t\ and t 2 and are given as 
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(3-1) 



(3-2) 
(3-3) 
(3-4) 



and tii(T) and t 2i (T) are regarded as the central values of the time of interactions and are 
defined by 
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(3-5) 



Because the explicit forms of ^(Td), t 2i (Tr>), a\ , af 2 . and At { are quite complicated, we 
just give quadratic forms of the trajectory function, and the explicit forms are given in 
Appendix [X] Here the central values of the times and time widths depend on the momenta 
of the external and internal particles. Therefore it happens in some situations that time 
widths become very large, although the spatial widths of external particles are finite. 
The amplitude Si is given by 



Si = N exp 
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(3-6) 



From the integrand in Eq. (J3-6J) . we find that the integrations over t\ and t 2 are separated 
when both time widths a\ and a\ are small enough compared to the time interval Tp. 
Afterwards, we assume that these conditions are satisfied. Then we can integrate Eq. (J3-6J) 
over ti and t 2 , and we obtain 



S; 



27r7V\/0f 1 . cr|. x exp 
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• (3-7) 



The quantity a/^L^L a PP ears as an overall factor in a consequence of the time integrations. 
The amplitudes become large when the time widths are large. This factor is derived only in 
a field theoretical treatment in which whole process is involved. 

In the second line of Eq. ([3-7)1 . the contribution from the trajectory function to the 
amplitudes becomes maximal at t^ and t 2i . Here, t^ is the central value of the production 
time of Jj, and t 2i is that of the detection time. The Gaussian integration with t^ and t 2i in 
Eq. f|3-6|) becomes negligible, unless the condition 



< t° ti (T D ) < t° 2i (T D ) < T D . 



(3-8) 



is satisfied. 

We see from Eq. (|3-7|) that the phase difference 2 \ up to 0{m 2 ) can be expressed in the 
form 



021 (T) 
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(3-9) 



m,i=0 , 



where t\, t 2 , AE® and AE 2 are the central times and energy differences with rrii = 0. 

The first term of Eq. (J3-9J) corresponds to the phase of the standard formula, and the 
second term results from the field theoretical treatment. The phase difference takes a form 
similar to that of the standard formula when AE®, AE 2 and AP are zero. In this case, we 
have the standard formula when t 2 — t\ can be regarded as the travel time of the neutrinos. 
However, t 2 is different from the final time To, and t\ is different from the initial time, and 
they are given by Eq. ()3-5|) . Thus, our formula is not exactly the same as the standard 
formula. 



The absolute square of the amplitude, IS") 2 , gives the transition probability from the 
initial state (|2-7|) . which is prepared at t = to the final state ()2-8j) . which is measured at 
t = Tp. In observations of solar neutrinos and atmospheric neutrinos, the detection time is 
not measured. In baseline neutrino experiments, the detection time is measured but with 
finite precision. Therefore we have to sum up IS" | 2 in a finite detection time. Here we make 
the assumption that the measurements at different times are independent phenomena. This 
assumption is consistent with the definition of wave packet states that they are defined with 
infinite precision and with the fact that the event rate is very small. From this assumption, 
we can integrate \S\ 2 over the detection time T D : 



\S\' 



dT D \S\ 2 . 



(3-10) 



The integration interval is different for different situations. It is from to oo for neutrinos 
from the sun or atmosphere, and from To — AT^/2 to To + ATp/2 for baseline neutrino 
oscillation experiments, where ATp is the resolution of the detection time. 

We assume that the Gaussian approximation is valid for the To integral. Then, the 
transition probability becomes 
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(3-iiJ 



where A t represents the energy and momentum functions and trajectory function, and o\ Di 
represents the width of the detection time, Tjj. These are given by 
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(3-12) 
(3-13) 
(3-14) 



The quantity 2 \ is the phase difference between 1^ and I 2 in Eq. (J3-9J) . The detection time 

is replaced by its central value, 

a 2 T° 4- rr 2 T° 
°T D l 1 D2 ^ °T D 2 1 Dl 



T, 
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a T D l + a T D 2 
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where T^ )i is the central value of the detection time in which only the mass eigenstate 7, 
appears as the intermediate state, and it is defined by 

dZ, 



= 0. (3-15) 

Td=T° Di 



dT D 

In the exponent of the second term in Eq. (|3-11|) . last two terms are characteristic terms 
in wave-packet treatment. One is called the "decoherence function", which represents an 
overlap in the detection time through two intermediate states : 

frpQ _ rpQ \2 

Decoherence function = — ^ ^ — -. (3-16) 

2(4 d1 + 4 d2 ) 

When the detection time difference, T^ 2 — T^ )1 , becomes larger than the detection time 



width, \ &t d i + ®t d 2i the oscillation disappears, because the coherence in the time direction 
is lost. The other is called the "phase function", which gives the condition that the detection 
time width must be smaller than the oscillation period : 

Phase function = -I*| (™22°> " 

2 Td \ dT r 



■D 

2 



l(Jr£s) (2tt-0 21 (O)) 2 . (3-17) 



Here T£f c is a period and is defined by 6>2i(Tg sc ) = 2tt. The extra coefficient 2n — A0(O) 
appears due to the field theoretical treatment. 



The energy function places a constraint on the energy differences AE^ and AE^, 



-i (° tli AEl) 2 - \ (^AE° 2i + At° u AEty 2 . (3-18) 

Energy conservation is satisfied when the energy function is zero. 

When the energy-momentum and the trajectory functions in one dimension are zero, the 
phase difference, Eq. (|3-9|) . becomes 

®2i(T° D ) = -^{Xb -X A + 4(T°)v b - t\(T» D )v A ). (3-19) 

Equation (|3-19|) agrees with the phase of the standard formula 9 )' 10 ' 1 if t°VB and t\v A are 
negligible. 

§4. The effect of a finite lifetime on the wave packets 

In this section, we study the transition amplitude and probability when a source particle 
has a finite lifetime. *' From these, the particle oscillation of the intermediate particle is 

*' This includes cases in which a particle is stable but the state looses the quantum mechanical coherence 
after a finite time. The r stands for the relaxation time in this case. 
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studied. 

The finite lifetime r is introduced for the particle A in the following manner. The field 
operator of the particle A, <£>a, contains r = 1/V as 

d 3 p 



$a(x) = I 2jjr exp (zp ■ x - i(E A - iT/2)t) a(p) A + h. 



(4-1) 



In consequence of this addition in Eq. (|2T2jl . a damping factor —T/2 is added to Ea(pa) 
in the previous section. In this section, we study the amplitude and probability in the case 
that both widths a t2 and a tl are small enough compared to the time interval To and the 
lifetime r satisfies r ^> a tl and a <2 . Then, as in the previous section, it has been found that 
the integrations over t\ and t 2 in the amplitudes are separated. 

If these conditions are not satisfied, the Gaussian approximation is invalid, and the 
integrations over t\ and t 2 cannot be separated. Then, the calculation would be similar to 
that for in plane waves. Below, we assume that these conditions are satisfied. Therefore the 
transition amplitude given by Eq. (J3-6J) are modified into the following form: 



where 



Si = 2nN\/(T$ 1 .(T% u x exp 



S = -F 1 F 2 sin29(S 1 -S 2 ), 
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(4-3) 



The new central times t^ and t 2i of the Gaussian functions of t\ and t 2 in Eq. (J4-3J) are 
obtained from old ones in Eq. (J3-5J) as follows 



t 2l '(T D ) = tUT D )-a 2l At° u ^. 



(4-4) 
(4-5) 



The transition probability is also calculated by assuming the Gaussian approximation for 
Td integration. It is given by 
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(4-6) 
(4-7) 
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where the coefficients Cj are the same as before. T^ )i is the central time of the detection 



time and is given by 



I Di~ 1 D i ~ -^ a T D i I q Td )• d ■•'■■) 



6>2i an d A^ are obtained by replacing £§ i; ^ an d Teh by £?j , ^ anc ^ ^d« m Eq. (|3-9|) and 
(J3-12J) . The Gaussian approximation for t\ and t 2 is useful and is valid when r ^> of., 
of is satisfied. This condition is satisfied for a lifetime r = 10~ 8 sec of the pion and 
o~ xL (L = A, B, C, D) of atomic size. 

We see from Eqs. (J3-11J) and (J4-6J1 that the exponential in the oscillation term becomes 
maximum at the given positions X, when the energy-momentum function and the trajectory 
function vanish. However the peak of the oscillation probability does not always coincide 
with that of the exponential, because of the coefficient and the lifetime. The lifetime and 
relaxation time reduce the magnitude of the oscillation probability depending on how long 
the particle A lives. If the constraints from energy-momentum conservation and the classical 
trajectories are weak, the oscillation probabilities become maximal at the position where the 
energy and momentum are not conserved. Then, as a result, the oscillation length changes 
from that of the standard formula. Such situations seem quite strange. However as we 
show below, these phenomena occur when the spatial widths of the external particles are 
extremely small, on the order of 1 fm. 

§5. Coherence conditions 

Here we examine necessary conditions for oscillation of an intermediate particles to take 
place, based on the transition probabilities Eq. (J3-11J) and (|4-6|) . Oscillation occurs when the 
interference term of two amplitudes becomes finite. Two amplitudes have peaks in different 
positions and decrease rapidly, so the interference term becomes finite only when the peak 
position overlaps within the widths. We find these conditions, coherent conditions, in the 
present section. There exist several previous works on this topic ? 3 > 7 ) 12 )> 13 ) but our results 
based on field theory are different from them. In our method, the coherence conditions are 
written in terms of measured quantities, like the positions, velocities and wave packet sizes 
for external particles. 

5.1. Energy function 
The factor 



exp 



2 1/ _, _ \ 2 



2 



(*««^?i) -^AEl + KAEl 



(5-1) 



13 



which also appears in the amplitude Si, yields the approximate energy conservation and 
imposes the constraint on Am\ x and a x to generate oscillations. The Gaussian function (J5-1J1 
has a width a\ of AEu and a\ of AE-n + -r^-AEu. Note that the second term in (J5TJ1 
contains AE U because of finite interaction time interval. We understand from these terms 
that when AE X 2 — AEu or AE22 — AE 2 \ becomes larger than (<Xt u ) -1 or (a t2i — Atu) -1 , the 
interference of I\ and I2 disappears, and the oscillation does not take place. From Eqs. (J2-24J1 
and (|2-25j) . AE 12 — AE U = ^km an d AE 2 2 — AE 2 ± = — 2 Tpl , this coherence condition for 
oscillations is reduced to 

Ami, 1 

21 < 



2E ~ \a tl \ 
Ami, 1 



'21 



< 1 3— r, (5-2) 

2E ~ \a t2 -Attl 

where E represents |k°| and a tl , a t2 and At\ are the values of Eq. ()3-2|) - (j3-4j) at m; = 0. 
The above coherence conditions (J5-2JI are expressed as 

2tt 1 1 . , 

j— <t—t, or — (5-3) 

where L osc = ^f- is the oscillation length in the standard oscillation formula, which is 
almost the same as that derived from (|3-19|) . Using the spatial widths a xi and a X2 of 
intermediate particles in the production and detection processes defined by 

&xi = vi\a tl \, 

&X2 =vi\a t2 - Ati\, (5-4) 

where vj is the velocity of the intermediate particle and is equal to 1, Eqs. (|5-3|) is written 

as 

&xi, &x 2 < -y^ (5-5) 

for standard oscillations to take place. When the above constraints are not satisfied, ordinary 
interference of I\ and I2 does not exist, and the oscillation disappears. Consequently, the 
transition probability becomes constant in the time interval Tjj. 

5.2. Trajectory function 

The trajectory functions in each amplitude, Si, 
1 



22 22 22 

QVT ^ a AC a BD T7i2Ai.O -0 rpQ \ °B°D r >2( + rpQ \ °A°Ctt2/ + 

za ^°BD Za AC 
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(5-6) 



also give necessary constraints for oscillation to take place. These forms are very complicated, 
and it is difficult to find expressions using parameters of external particles. For this reason, 
we express constraints using the center times and the center positions. These conditions are 
given as 



(5-7) 
(5-8) 
(5-9) 




° 2 A°l 



From Eqs. (|2Tfij) to (|2T8j) . the above constraints are rewritten as 



*x° - 64 - k(Sf 2 - 6%) + 4#k(t°(T° J - tUn J) 



< 



a z 



2E 2 

\(v B -v D )5t° 2 + v D 5T°\ < 

Wc - v a||^i| < 



a AC a BD 



a 



BD 



2 2 ' 

°b°d 



a 



AC 



°A°c 



(5-10) 
(5-11) 
(5-12) 



Throughout this section, we write the center times without a prime for the case in which 
the particle A has a finite lifetime. 

The right-hand side of the first constraint, Eq. (J5T0J1 is the spatial width of the inter- 
mediate particle, and those of the second and third constraints, Eq. (J5-11J) and (J5-12J1 . are 
the sums of the spatial widths of the external particles. On the left-hand side of the first 
condition, ^l 1 k(tj] — t?) is rewritten as follows 



> -^travel f 

A— k, 



(5-13) 



where L trave i = t\ — t®, and A is the de Broglie wavelength for the intermediate particles. 
Unless ^travel is much larger than L osc , Eq. ()5-13|) is of order A or less, and the condition 
(|5T0j) is satisfied. 

5.3. Decoherence function 

The decoherence function Eq. (I3T6J) appearing in the probabilities (J3-11JI and ()4-6|) con- 
strains the difference between the central times of detection time as 

5T D ° < 2a T . (5-14) 
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Here, &t d is the value of Eq. (|3T4|) at rrii = 0. This constraint gives the "decoherence time" 
as a function of the detector position, particles velocities and wave packet sizes. Note that in 
our approach, oscillating particles appear as intermediate, states and only external particles 
are observed. Therefore the decoherence condition is given for the detection times of the 
scattering particle, not for the flight distance of intermediate particles. From (|5T4jl . when 
STd° is larger than 2ot d -, coherence is lost and no oscillation is seen. 

5.4. Phase function 

The phase function Eq. (J3T7JI in the oscillation probability gives a constraint on oscilla- 
tion period, Tg sc : 

2(2n-O 21 (0))a TD <T° D sc . (5-15) 

This relation implies that when the oscillation period is smaller than the width of the detec- 
tion time, the oscillation disappears. 

5.5. Lifetime 

The lifetime effect from a source particle A is seen explicitly as exp(— r/2tf) in the 
absolute square of the amplitude (J4-6J) . From the right-hand side of (J4-7J) . this term gives 
constraint 

(ft? < r. (5-16) 

From this relation, to maintain coherence, the difference in the production time for interme- 
diate particles must be smaller than the lifetime of the source particle. 

§6. The numerical results of transition probabilities 

In this section, we give the results for the numerical calculations of the oscillation proba- 
bilities. The particle C, which corresponds to a muon in pion decay, usually is not detected 
in most experiments and observations. From this fact, the oscillation probabilities that we 
actually measure are the sums of probabilities over p^. Therefore, we consider the following 
probabilities instead**) of those in Eq. ()3T1|) and Eq. (|4-6|) : 

P(X B ) = N novm [dp° c \S\ 2 . (6-1) 



**> Actually, \S\ 2 is constant with Xc in finite macroscopic range. Therefore this value is equivalent to a 
probability integrated over coordinate 

Jdp° c JdX c \S\ 2 , 

in which the orthogonality of states with different values of pc° is satisfied. 
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Here, \S\ 2 is given in Eq. (J3-11J) or (|4-fij) . and iV norm is a normalization factor. We investigate 
the following three situations of different parameters: 

1. Intermediate particles are produced by the decay of the particle A in flight. The average 
momentum of each external particle is taken so that the average momentum of the 
intermediate particle is about 430 MeV. This case mimics long base line experiments. 

2. The intermediate particles are produced by the decay of the particle A at rest. The 
mass of the particle A is chosen as about 140 MeV, which gives the intermediate 
particles' momentum as 30 MeV. 

3. The intermediate particles are produced by the decay of the heavy particle A in flight. 
The average momenta of external particles are taken so that the intermediate particles 
momentum becomes 0.4 MeV. 

The first and the second cases correspond to "decay in flight" (DIF) and "decay at rest" 
(DAR) neutrinos in neutrino oscillation experiments using pion sources. The third case 
correspond to "solar neutrinos from 7 Be decay" whose energies are MeV (low energy or LE). 

In baseline experiments, the momenta of source particles are focused in the direction of 
the detector and the momenta of produced neutrinos and accompanying charged leptons are 
in almost the same direction. In solar neutrino observations, the velocities of source and 
accompanying particles are slow, because their momenta are much lower than their masses. 
Therefore the one-dimensional approximation is valid in the above three cases. For these 
reasons, we perform p^ integral in one dimension. 

For the numerical calculations, we use the following parameters: m x = 0.1 eV and 
m 2 = 0.2 eV and 9 = ~ in all cases. These are the same values in (I). The wave packet 
sizes, o"l (L = A,B,C,D), are also taken to be the same. The values of other parameters 
are given in Table |U 

Tc in case 1 and case 2 and Xq in case 3 are not shown in Table U] We set these 
parameters by hand as a function of p^,, because there is no way to determine both Xq and 
Tc in one dimension. The concrete forms of Tcip%) and Xc(pc) are given in the following 
subsections. 

6.1. Case 1: Decay inflight 

In the first case, the particle / is produced by the decay of the particle A in flight 
accompanying C, and the particle D appears through the interactions between the particle 
B at rest and /. 

In baseline experiments, the particle C is considered to be stopped at a beam dump. 
Therefore we set Xq to a constant value. But the detection time of C is unknown. Therefore 
we study the probability of a certain time T c , which is a function of Pq. T c should be such 
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lifetimes or relaxation times, r are in seconds. (L = A, B, C, D) 



that the time order for the central times is given by 



o < und < UTid < T" Di 



tl < Tc, 



where t\^ t^ and T^ )i are functions of T c . One choice satisfying this condition is 



Tc(p°c 



X, 



c 



C) 



vc(p°c) 



x 0.999. 



(6-2) 
(6-3) 



(6-4) 



Using Eqs. ()6-4|) and (j6-l|) . the oscillation probabilities are calculated numerically. 

The oscillation probabilities with an infinite lifetime and a finite lifetime are shown in 
FigsHJandini respectively. In Figs. El and HH the oscillation length becomes longer than that 
of the standard formula when the wave packet sizes are smaller than 3.0 x 10~ 15 m, and the 
amplitude of the oscillation probability becomes smaller than 1 as X B becomes larger or as 
the wave packet sizes become smaller. 

In Figs. 0] (a), (b), (c) and (d), we compare the oscillation probabilities of the infinite 
lifetime with those of the finite lifetime at a x = 2.0 x 10 -15 m for (a) and (b) and a x = 
1.0 x 10~ 15 m for (c) and (d). From Fig. HJ it is seen that the oscillation length is longer and 
the amplitude of oscillation is smaller when the lifetime is infinite. 
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1.0*10" 15 m 
2.0*10" 15 m 
3.0*10" 15 m 
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X B (km) 



Fig. 2. The DIF oscillation probability with infinite lifetime. The solid curve represents the stan- 
dard formula (S.F.), and the dashed, dotted and dashed-dotted curves correspond to the cases 
in which the wave packet sizes are 1.0 x 10~ 15 m, 2.0 x 10 -15 m and 3.0 x 10~ 15 m, respec- 
tively.The horizontal axis is the position of B, Xb (km). 



S.F. 




1.0*10" 15 m 
2.0*10" 15 m 
3.0*10" 15 m 



120 



X B (km) 



Fig. 3. The DIF oscillation probability with finite lifetime. The solid curve represents the standard 
formula (S.F.), and the dashed, dotted and dashed-dotted curves correspond to the cases the 
wave packet sizes are 1.0 x 10~ 15 m, 2.0 x 10~ 15 m and 3.0 x 10 -15 m, respectively. The 
horizontal axis is the position of B, Xb (km). 

The increase of oscillation length is caused by the increase of the time widths, <r 4l . and 
OTii in the coefficients in Eq. ()3T1|) . In the following, we clarify the reason that the oscillation 
length increases with the time widths. 

In Fig. the p° c dependences of the time widths are shown. From Figs. El (a) and (b), 
it is seen that a tli and a r . grow with p° c . From Eq. ()2-18|) . it is seen that H loses its ti 



19 



_Q 

CO 
_Q 
O 



q 
jo 
o 

CO 

o 



-O 

CO 
_Q 
O 



g 
jo 
o 

CO 

o 



1 




1 


A ' 


1 


'A 






1 


A '* 


1 


■ 


1 


1 




0.9 










A (a) ■ 






0.9 










(b) - 




0.8 










7 \ j 




>. 


0.8 






J u 








0.7 














s 


0.7 






\ ll Ik 
i 7 \ '/1 








0.6 
0.5 
0.4 
0.3 
0.2 






» 

A 
A 


3 

1/ 
■ 


M '/ 
l\ J - 





a! 

.0 


CO 




0.6 

0.5 
0.4 
0.3 
0.2 






1 1 1 a \ 
\ I \ i \ 

Si 11 Af 

5 1 \ 2 A/ 


1/ 

A3 


\ A A 

Ay v 

vy - 




0.1 


( 






A 

\\ 
\ 










0.1 



c 














) 20 


40 


60 


80 


100 120 


) 50 


100 


150 200 


250 


300 350 








X B (km) 














X B (km) 






1 

0.9 
0.8 














>. 


1 

0.9 
0.8 


















// \ 
7 \ 




A " : 












(d)- 




0.7 










A 




s 


0.7 














0.6 










'/ \ 




O 


0.6 






/\ 








0.5 
0.4 
0.3 
0.2 






A 
A 

n 
1 
1 


\a 


ft j( 




CL 

.0 


CO 




0.5 
0.4 
0.3 
0.2 




f \ i 
1 " 


(V^v 














0.1 








\j 








0.1 


' 













( 


/ 




, 


, 


„ 







C 




, 


. 








) 20 


40 


60 


80 


100 1: 


>0 


) 50 


100 


150 200 


250 


300 350 








X B (km) 
















X B (km) 









Fig. 4. The oscillation probabilities at a x = 2.0 x 10 -15 m in (a) and (b), and those at a x = 
1.0 x 10~ 15 m in (c) and (d). The solid curves represent the oscillation probabilities with 
infinite lifetime and the dashed curves represent those with finite lifetime. 



dependence when vc equals va- As we mentioned in the last part of Section 4, this growth 
of the time widths shifts the peak of the oscillation probability upward from that of the 
Gaussian in Eq. (JHTlj) . 

FigurelHl(a) shows the Gaussian function exp(2v4j) in the infinite lifetime case, and Fig. EI 
(b) shows the absolute square of total amplitude Eq. (|3T1|) . in which the maximum values 
are normalized to unity, at a = 1.0, 2.0, 3.0, 4.0 x 10~ 15 m, respectively. From Fig. H (a), 
it is seen that the peaks of the Gaussian function are independent of the wave packet sizes. 
From Fig. © (b), it is seen that the peaks of the oscillation probability depend on the wave 
packet sizes. For large a x (> 4.0 x 10~ 15 m), the peaks of the oscillation probability are 
the same. But for small a x (< 4.0 x 10~ 15 m), the peaks of the oscillation probability are 
different from those of the Gaussian functions. This happens because the time widths grow 
with Pq in Fig. EJ and the width of the Gaussian function becomes large in Fig. El (a). Then, 
the main contribution to the oscillation phase, Eq. (J3-9J1 in Eq. ()6Tj) . comes from the higher 
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Fig. 5. (a): The p° c dependence of a% x . The solid curve corresponds to <r tl at o x = 1.0 x 10~ 15 m, 
and the dashed curve corresponds to a x = 3.0 x 10 -15 m. (b) :The p° c dependence of <7 i2 and 
ot- The solid and dotted curves represent ct 42 at a x = 1.0 x 10 -15 m and a x = 3.0 x 10~ 15 m, 
and the dashed and dashed-dotted curves are &t at a x = 1.0 x 10~ 15 m and a x = 3.0 x 10~ 15 
m. 
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Fig. 6. (a): The p° c dependence of the Gaussian function, (b): The p° c dependence of the first 
term in Eq. (|3T1|) . The solid curve represents a = 1.0 x 10 -15 m, and the dashed, dotted and 
dashed-dotted curves represent 2.0 x 10 -15 , 3.0 x 10 -15 and 4.0 x 10~ 15 m, respectively, in both 
graphs. 

p° c region. 

Figure [7| displays the p Q c dependence of the ratio of the phase difference Eq. ()3-9|) to that 
of the standard formula for Xb = 20, 60 and 90 km : 



[Am 2 



Phase Ratio = 21 {f %)/ —^X B 

\ Ah 



(6-5) 



Here Xb is regarded as being the same as the travel distance of intermediate particles, 
and the energy E in the phase of the standard formula is determined by AP = and 
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Fig. 7. The p° c dependence of the ratio of oscillation phases. The solid, dashed and dotted curves 
represent the ratios at Xb = 20, 60, 90 km, respectively. 



AE® ~ AE\ — 0. From Fig. [7[ it is seen that the ratio is almost independent of Xb but 
depends on p° c . As is shown in Eq. (I3-19J1 . the phase ratio should become 1 at the solution 
AP = and AE\ = AE% = 0, p° c = 570 MeV. When a x is larger than 3.0 x 10~ 15 m, the 
main contribution to the integral in Eq. (|6T|) comes from the region where the energy and 
momentum are conserved. Equation ()6T|) is almost the same as the standard formula. But 
from Fig. H3 when a x is smaller than 3.0 x 10 -15 m, the peak of the Gaussian part in the 
integrand moves to higher p° c region than that of Fig. |U] (a), and the main contribution to 
Eq. f|6-l|) comes from the higher p° c region. Then, the oscillation length becomes longer than 
that of the standard formula. 

In the case that the particle A has a finite lifetime, the main contribution to the integral 
comes from the lower p° c region, because of the presence of the lifetime in Eq. (J4-7J) . 

Figure |H1 (a) displays the p° c dependence of exp(— rt\.), and (b) displays the p° c de- 
pendence of the absolute square of the amplitude (|4-fij) . in which the maximum values are 
normalized to unity, at o = 1.0, 2.0, 3.0, 4.0 x 10~ 15 m. The quantity exp(— -TY?.) is almost 
independent of X B and a x . From Fig. El (b), the peaks of the oscillation probability Eq. (J4-6J) 
shift to the lower p° c region as the wave packet sizes become small. 

The phase ratio Eq. (J6-5J) with finite lifetime is almost the same as that with infinite 
lifetime, because the lifetime r is much larger than the time widths. Therefore, from Fig. 
the oscillation length is longer than that of the standard formula when the wave packet sizes 
are smaller than 2.0 x 10~ 15 m and are slightly smaller than that of the standard formula 



when a x is larger than 3.0 x 10 m. 

Figure M displays the p° c dependence of the absolute square of the flavor changing ampli- 
tudes around X B = 37 km. It is seen that the peaks of the absolute square of the amplitudes 
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Fig. 8. (a): exp(— -Ti^) in the exponent of Eq. (|4-6j) . (b): The first term in Eq. (|4-6j) . The solid and 
dashed curves correspond to 1.0 x 10 -15 m and 2.0 x 10 -15 m, and the dotted and dashed-dotted 
curves correspond to 3.0 x 10~ 15 m and 4.0 x 10~ 15 m, respectively. 



change and the shape of each amplitude is deformed with distance. The Gaussian shapes 
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Fig. 9. p® c dependence of the absolute square of the flavor changing amplitude at a x = 1.0 x 10 -15 
m. The vertical axis is normalized by the value of \S\ 2 at Xb = 18 km. The solid curve 
represents the flavor conserving process at Xb = 35 km, and the dashed and dotted curves 
represent the flavor changing processes at Xb = 37, 39 km. 
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of the amplitudes are deformed by the cosine in the oscillation term. The behavior of the 
oscillation depends on Xb as well as p° c , since the phase is proportional to X^/|k |. Then, 
as a result, the change and deformation given in Fig. |U] occur. This effect is seen clearly 
around the oscillation minimum. 
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6.2. Case 2 : Decay at rest 

In the DAR, the intermediate particles I and particle C are produced by the decay of 
particle A at rest. The particle / is scattered by the particle B at rest, and then the particle 
D appears. 

For the same reason as in case 1, Xc is set to —5.0 m and Tc is given as 

X c 



Tc(p°c 



c) 



voipc 



x 2.0. 



(6-6) 



C) 



Because the source particle A is at rest, the numerical factor in Tc{pc) must be greater than 
1 to satisfy the conditions (J6-2JI and (|6-3|) . Here we take this factor to be 2. 

The oscillation probabilities with finite lifetime are shown in Fig. EB It is seen again that 
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Fig. 10. The DAR oscillation probability with infinite lifetime: The solid curve represents the 
standard formula (S.F.), and the dashed, dotted and dashed-dotted curves correspond to the 
cases in which the wave packet sizes are 2.0 x 10 -14 m, 3.0 x 10 -14 m and 4.0 x 10 -14 m, 
respectively. The horizontal axis is the position of B, Xb (km). 

the amplitude of the oscillation becomes smaller and the period of the oscillation probability 
becomes longer than that of the standard formula as the wave packet sizes become smaller 
than 4.0 x 10~ 14 m. 

In this case, t^ is almost zero, since the particle A is at rest. Because of this, there are 
no significant differences between the source particle with a finite lifetime and an infinite 
lifetime. 

6.3. Case 3 : Low energy 

In the last case, the intermediate particles are produced by the decay of the heavy particle 
A in flight. The central value of its momentum is about 1.3 MeV, which is much lower than 
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in case 1. The particles A and C have larger masses and smaller momenta than in the other 
two cases. 

In contrast to the above two cases, this case corresponds to the solar neutrinos from 
7 Be decays. The change of quantum mechanical states by scatterings is considered to be 
equivalent to detection or observation. Therefore the detection time of C is taken to the 
relaxation time, which is assumed to be the same value of particle A, and Xc is given as a 
function of p c - The quantities Tc and Xc are 



T c 



X c (p c ) = 1.5 x T c v c (Pc), 
where r is the relaxation time of the particle A, and its value is set to 10 -12 sec 



(6-7) 
(6-8) 



In Fig. ^TJ the oscillation probability with finite lifetime is shown. In this case, the 
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Fig. 11. The LOW oscillation probability with lifetime: The solid curve represents the standard 
formula (S.F.), and the dashed, dotted and dashed-dotted curves correspond to the case in 
which the wave packet sizes are 1.0 x 10~ 13 m, 2.0 x 10~ 13 m and 3.0 x 10~ 13 m, respectively. 
The horizontal axis is the position of B, Xb (m). 

oscillation length becomes longer and the oscillation amplitude becomes smaller again. The 
wave packet effect becomes observable if the wave packet sizes are of order 10~ 13 m. A wave 
packet size of order 10~ 13 m is larger than those for the DIF and DAR cases. This is because 
the central values of the momenta are lower than in the DIF and DAR cases. 

§7. Discussion and conclusion 

In this paper, we have studied particle oscillations in the intermediate state of transition 
amplitudes of wave packets based on a simple scalar model. A source particle, a target 
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particle, and two scattered particles are the wave packets. In this situation, wave functions 
do not spread infinitely, but are localized within a finite width a x , and interference occurs 
only among intermediate particles that overlap spatially. This interference disappears and 
the oscillation probability deviates from those of the standard formula in parameter regions 
where the intermediate particles are separated spatially. 

We computed the total oscillation probability and the phase factor of the amplitude 
numerically. We found that the oscillation probability agrees with that of the standard 
formula if the wave packet sizes are of semi-macroscopic values, and that the oscillation 
probability deviates from those of the standard formula in extreme parameter regions. This 
occurs when the wave packet sizes are of the order of 10~ 13 m or smaller. In this region, 
the wave packet size effects become visible. The oscillation amplitude becomes smaller and 
the oscillation period becomes larger than those of the standard formula. We hope that 
this region may be realized experimentally and that the modified formula found here will be 
tested in the future. 

Although our results were obtained on the basis of a simple scalar model, we hope that 
for precision measurements of the neutrino parameters such as masses, the MNS matrix, and 
others, our considerations will be valuable. 
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Appendix A 

The Center Times and the Time Widths 

In this appendix, we give the explicit forms of the time widths (af u ,af 2 .,<f^ Di ) and the 
center times (t^, ^ij^dJ- We omit the mass index % for simplicity. However, we note that 
one can easily obtain the center times and the time widths for specific mass eigenstates by 
replacing vj with z/j. 

The time widths and the center times are written in terms of the time derivatives of the 
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classical trajectory. The classical trajectory given in Eq. (|2-15jl is 



Zi(ti, t 2 ,T D ) 
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We use following abbreviated expressions for the time derivatives of the classical trajectory 
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Here, Fj, Gj and Hj [i = 0, 1, 2, T) are the coefficient vectors of ti, t 2 and T, and are given 



as 
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F 1 = - ^ V ^ VC +v„ F^ ^ + ^ -v,, F T = -^, (A-15) 

^AC a BD °BD 

G = X D - X B , G 2 = -v B + vd, Gt = -v D , (A-16) 

H = X c - X A - T c v c , K 1 = -y A + v c . (A-17) 

Using these functions, the time widths and the center times are written as follows : 



— = ^n, (A-18) 
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_ = _ Ztt - -(a 2 ti + At\)Z[ T - -a t2 Af 1 Z 1T Z 2T - -af 2 Z 2T , (A-21) 



*i(r D ) = -^(4 + ^?Vio - ^AfiZn, (A-22) 

*2^) = —^At^o - \o\Z 2 ^ (A-23) 

T° = -^4Z T0 - (A-24) 

Appendix B 

— Measurement of Transition Probability 



Following the standard interpretation of measurement in quantum mechanics, the square 
of the absolute value of the amplitude, \S\ 2 , gives the transition probability from the initial 
state (J2-7J) . which is prepared at t — 0, to the final state Eq. (J2-8J) . which is defined at t — Tr>. 
For the probability interpretation to make sense, each observation at the final state should be 
made independently and exclusively. This is, when one value is observed for an observation, 
the other value should not be observed. The state of one value is different from the state of 
a different value. Conversely, the final state should be different if the corresponding values 
are different. 

Because the orthogonality of states described by the wave packet is peculiar, the total 
probability should be defined in a manner that is consistent with experiments. This problem 
has been solved in usual scatterings, where the wave packet effects are negligibly small. 
The detector used in experiments has a finite macroscopic size, and the total rate observed 
in a macroscopic detector is computed with a continuous momentum. The total observed 
probability is obtained by integrating the square of the absolute value of the amplitude with 
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a continuous momentum. Therefore this problem has been solved in the case of normal 
scattering. 

In neutrino scattering, the energy scale is very small, and the event rate is also very 
small. Events occur so infrequently that the detector is in different quantum states that 
are orthogonal to each other. Also, the detector has a macroscopic size. Consequently, the 
total rate observed in a macroscopic detector within a finite detection time is computed by 
integrating the probability with the central value of the momentum and the time. 
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